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Photosynthetic light harvesting provides a nat- 
ural blueprint for bioengineered and biomimetic 
solar energy and light detection technologies. 
Recent evidence |ll4ll| suggests some individ- 
ual light harvesting protein complexes (LHCs) 
and LHC subunits efficiently transfer excitons 
towards chemical reaction centers (RCs) via an 
interplay between excitonic quantum coherence, 
resonant protein vibrations, and thermal decoher- 
ence. The role of coherence in vivo is unclear 
however, where excitons are transferred through 
multi-LHC/RC aggregates over distances typi- 
cally large compared with intra-LHC scales |12i - 
Here we assess the possibility of long-range 
coherent transfer in a simple chromophore net- 
work with disordered site and transfer coupling 
energies. Through renormalization we find that, 
surprisingly, decoherence is diminished at larger 
scales, and long-range coherence is facilitated by 
chromophoric clustering. Conversely, static dis- 
order in the site energies grows with length scale, 
forcing localization. Our results suggest sustained 
coherent exciton transfer may be possible over 
distances large compared with nearest-neighbour 
(n-n) chromophore separations, at physiologi- 
cal temperatures, in a clustered network with 
small static disorder. This may support findings 
suggesting long-range coherence in algal chloro- 
plasts [2], and provides a framework for engi- 
neering large chromophore or quantum dot high- 
temperature exciton transfer networks. 

Individual LHCs exploit a balance of quantum coher- 
ence and thermal decoherence to improve their exciton 
transfer efficiency. At intra-LHC scales, vibrational reso- 
nances in the surrounding protein help to preserve exci- 
tonic coherence when electronic transfer couplings and 
protein reorganization energies converge [H, 0, 0, EH- 
Thermal dephasing noise cooperates with this coherence 
to transfer excitons through the network by suppress- 
ing transfer-inhibiting interference effects and exploit- 
ing spectral broadening of site energies to overcome An- 
derson localization; the result is transfer efficiency ex- 
ceeding that possible thr oug h purely coherent or purely 



incoherent mechanisms jlfj, An interesting ques- 



tion is whether these effects can improve physiologically 
relevant, long-range transfer efficiency through multi- 
LHC/RC networks and, if so, whether this enhancement 
is due only to accumulated intra-LHC enhancements or 



some more subtle long-range mechanism. 

An arbitrary chromophore network with disordered 
coupling topology may be considered as a hierarchy of 
nested clusters, where a cluster is defined as a group of 
sites coupled more strongly to each other than to sites 
outside the cluster, and the degree of clustering is quan- 
tified by the ratio of intracluster to intercluster coupling 
energies. At the smallest scale are clusters of strongly 
coupled chromophores, such as monomeric LHC subunits 
(Figure[T]f). At the next scale, each chromophore cluster 
forms an effective site, which itself clusters with other 
similar sites; monomeric LHC subunits cluster to form 
whole LHCs. The hierarchy may proceed further, to 
LHC aggregates and/or LHC-RC supercomplexes, to ag- 
gregates of supercomplexes, which in some cases finally 
cluster in 3D membrane stacks called grana [l6[ . The ef- 
ficiency of transfer to RCs depends strongly on structure 
at multiple scales within this hierarchy. Within some 
species, organization of LHC aggregates can facilitate 
functional switching between a high transfer efficiency 
state and a so-called nonphotochemical quenching state, 
which dissipates excess excitons as heat [l7j]. However, 
mechanisms of this switch remain controversial and their 
relationship to fundamental exciton transfer mechanisms 
is yet to be fully understood. Chromophoric clustering 
hierarchies also exist in some artificial li ght harvesting 
materials such as organic semiconductors |l8l . [l9| . 



Dimer 

Initially we study the crossover from coherent (tun- 
nelling) to incoherent (hopping) energy transfer in a 
pair (dimer) of undriven chromophores with separation 
d = |d| = |r( 2 ' — rW|. This dimer is our basic clustering 
unit, loosely analogous to a LHC subunit. We assume 
for the complete Hamiltonian [20j | . 



H = H- 



H 1 



H 



SB 



where H s is a spin-boson type, two-site transfer Hamil- 
tonian with site energies E^ 1,2 \ static disorder e = 
E^ 1 ' — E( 2 \ and transfer coupling A. The bath Hamil- 
tonian H B = fi-^q W ( c l) a ( c l)^ a (pi) is the sum over har- 
monic vibrational modes with wave vectors q and fre- 
quencies w(q). The system-bath interaction is H SB = 

*EqG (12) (q) (a(q) - a(-q) f ), where G( 12 )(q) is an ef- 
fective exciton-phonon coupling operator for the dimer 
(see Supplementary Information). 
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FIG. 1: Structural heirarchy of the light harvesting 
machinery from higher plants, a) Transmission electron 
micrograph of a sugarcane chloroplast, fixed, sectioned and 
stained. Densely stacked grana regions of the thylakoid mem- 
brane are visible. Image courtesy of Prof. Robert Birch, 
used with permission, b) Pseudo-3D cross-section of a granal 
membrane stack showing embedded aggregates of photosys- 
tem II-light harvesting complex II (PSII-LHCII) supercom- 
plexes interspersed with other complexes, c) 3D schematic 
of membrane-embedded semicrystalline aggregate of PSII- 
LHCII supercomplexes. d) PSII-LHCII supercomplex con- 
taining LHCs (green, brown) coupled to RC-containing core 
complexes (blue), e) Trimeric LHCII. f) Monomeric sub- 
unit of LHCII, with crystal structure overlaid. This complex 
binds 14 chlorophyll (8 chlorophyll-a, 6 chlorophyll-b) and 4 
carotenoid chromophores. g) Chlorophyll-a chromophore. 



We derive the system-bath coupling, 5 ( - 12 ^(q), spec- 
tral density and a Markovian master equation for pi, 
the combined density matrix of the site populations 
in an interaction picture with respect to H s + H B 
(see Supplementary Information). The spectral den- 
sity, J(uj) = 27r^ q <5(u; — tj(q))|g( 12 )(q)| 2 , which quanti- 
fies the strength of the system-bath coupling, is 



J(u) =B[uj i sin — 

a \ v 



(1) 



where B = HD 2 /(4tt/j,v 5 ) with D the exciton-phonon de- 
formation potential, fi the mass density and v the speed 
of sound. In the ordered case (e = 0), the site populations 
of an initially localized exciton show coherent tunnelling 
between sites at frequency 2A, damped at temperature- 
dependent decoherence rate F(2A). We call the under- 
damped regime, where quantum tunelling dominates, 
'coherent', and the critically- and over-damped regime, 
where thermal dephasing and hopping dominate, 'inco- 
herent'. The two regimes are defined as 

Coherent: F(2A) < 2A, Incoherent: F(2A) > 2A. (2) 



Clustered network and renormalization 

Consider a network of coupled dimers extending to 
length scales large compared with the intradimer sep- 
aration, analogous to an aggregate of elementary chro- 
mophore clusters such as LHC subunits (Figure [2]) . Our 
treatment is extensible to an arbitrary network geometry 
but we assume a one-dimensional arrangement, symmet- 
ric about the origin, with uniform intradimer separations 
do and interdimer (centre-centre) separations d\ such 
that d\ > 2do. This forms the 0th level of the structural 
hierarchy, which we now analyze through renormaliza- 
tion. We relax the assumption prevalent in earlier studies 
, l2ll I22I ] , that there is some length scale at which 
tunnelling becomes unimportant. Instead, we consider 
the dynamics at higher hierarchy levels as perturbations 
to the dynamics at lower levels [23j , allowing for the pos- 
sibility of coherent transfer at all scales. To assess the dy- 
namics at each level, we formulate an iterable renormal- 
ization procedure by which a cluster of n sites at the fcth 
level is mapped to a single effective site at the (fc + l)th 
level, located at r^+i = i J2]=i r k^ (See Supplemen- 
tary Information) . This is achieved via a coarse-graining 
operator which maps the n single-excitation electronic 
eigenstates (excitons) of the cluster to the unique single- 
excitation cigenstatc of the effective site. 

Figure [2] shows the renormalization flow for hierarchy 
levels 0-2. The flow is characterized by transformations of 
the site, transfer coupling and system-bath coupling en- 
ergies between levels. At the fcth level, the excitation en- 
ergy of a site is approximately the average of its Oth-level 
intracluster site energies and the transfer coupling be- 
tween sites is approximately the average of the Oth-level 
intercluster transfer couplings (see Supplementary Equa- 
tions (11) and (13). The system-bath coupling at the fcth 
level is given by the exciton-phonon coupling strength, 
M(q), modulated by a product of periodic terms depen- 
dent on the network geometry according to 



(12) 

9 k 



(q) = M(q) sin 



ifc • q 



fc-i 

n< 



d, q 



(3) 



The spectral density for an effective dimer at the fcth (fc > 
1) hierarchy level, </fc(w), is found by letting g^ 12 \q) — > 

( 12) 

g y k ' (q) in the spectral density and substituting equation 
©. This yields 



J fc (w) = BW k {u)u 3 



(4) 



where Wfc(w) is a modulating function that depends 
on the network geometry and characterizes the changes 
in the bath's coupling to the network at different lev- 
els (see Supplementary Equation (10)). The coefficients 
Ai(oS) = udi/v, with di = 2( l ~ 1 )di (i > 1) the intradimer 
separation at the zth level. The effect of the modulat- 
ing function Wk{uj) is to filter the system-bath interac- 
tion such that the phonon wavelengths which strongly 
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FIG. 2: Renormalization flow of dimeric chromophore network hierarchy. A tetramer at the Oth level is coarse-grained 
to an effective dimer at the 1st level via the operator . This dimer is paired with a similar effective dimer and the coarse- 

(i) 

graining repeated to obtain an effective dimer at the 2nd level, and so on. The site energies E k , transfer coupling energies 
A^' and system-bath coupling energies p^ 12 ' are transformed between levels. At each level, analogous structures from the 
higher plant photosynthetic structural hierarchy are presented in parallel with the dimeric network hierarchy: 0) Chlorophyll-a, 
1) LHCII monomer, 2) LHCII trimer. 



affect the network at a given level are proportional to 
the length scales characteristic of that level. At suffi- 
ciently high levels, the transfer coupling is diminished 
such that the small- frequency limit u <C v/d k holds and 
Jfc(w) ps Bd? k ui 5 1 (6v 2 ), from which it follows that 



dk 



k+l 



Jfe(w). 



(5) 



In this limit, given our chosen network geometry, 
Jfc+i (w) ps 4J fc (cj) for k > 1. 
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Scale-based crossover 



We assess the transfer dynamics at the fcth level in an 
ordered network (eg = 0) by calculating the ratio of 
the decoherence rate to the transfer coupling frequency, 
F fc (2A fc )/(2A fc ) where F k (2A k ) = J k (2A k )(N(2A k )+±), 
where N(u) = ( e hul / kBT — is the phonon occupa- 
tion of the bath at temperature T, and compare this with 
definitions in equation ([2]) . Figure [3] plots this ratio for 
the first six levels in ordered networks with varying de- 
grees of clustering (di/dn > 2). The parameter set used 
to generate figure El chosen for biological relevance (see 
Supplementary Information), supports coherent transfer 
at hierarchy levels 0-3 for all degrees of clustering, since 
Ffc(2Afc)/(2Afe) < 1. However, at levels 4 and 5 transfer 
is coherent only in more clustered networks, suggesting 
that clustering facilitates coherent transfer rather than 
inhibiting it, as may be assumed. Although decoherence 
dominates at intermediate hierarchy levels with low de- 
grees of clustering, at high levels (where 2A& <C v/dk), 
F fc+1 (2A fc+1 )/(2A fc+1 ) ~ 4 (Afc + i/Afc) 3 i 7 ? i ;(2Afe)/(2Afc), 



FIG. 3: Scale-based crossover simulation results. 

Shown for first six hierarchy levels (curves labelled) at 
T=293K in ordered (e^ = 0) networks with various degrees 
of chromophoric clustering (di/do). Dashed line separates co- 
herent (below) and incoherent (above) dynamical regimes. 



and therefore at sufficiently high levels, Ffe(2Afe)/(2Afc) 
vanishes. This indicates that in a large, ordered net- 
work, coherent transfer is supported at the largest length 
scales accessible within the excitation lifetime, provided 
it is also supported at all smaller scales. The latter re- 
quirement is met when the transfer coupling occupies the 
small frequency regime at all length scales. 

We interpret these results in the following way, assum- 
ing rapid transfer with sustained coherence may be bene- 
ficial to photosynthetic efficiency [l - 11 1 . In a biologically 
relevant parameter space of our network, rapid coher- 
ent oscillations dominate at small scales but are quickly 
damped. While coherence is more sustained at larger 
scales, the transfer rate is disadvantageously slow. Chro- 
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mophoric clustering balances these two extremes at in- 
termediate scales by restricting the bath modes able to 
couple to the network, thereby reducing the decoherence 
rate relative to the transfer rate. Consequently, coherent 
transfer is supported at larger length scales than in an 
equivalent unclustered network. 



Static disorder forces localization 



We now introduce intradimeric static disorder at the 



Oth hierarchy level, e, 



E 



(i) 



4 j) (* odd, j 



i + 1), sampled from a Gaussian distribution with vari- 
ance a 2 (eQ^). Since Gaussian variances add directly, the 
variance and mean of the static disorder at the fcth hier- 
archy level respectively are 



2 k- 



a e 



(4 IJ) ) 



The mean intradimer static disorder increases with k. At 
some level, ftA, this disorder will exceed the tunnelling 
rate (which decreases with k), i.e. e^f' >• , so 

the excitons become localized on a single effective site. 
This is a manifestation of Anderson localization. Previ- 
ous studies , |Tl| have found that at the Oth hierarchy 
level, spectral broadening due to thermal dynamic disor- 
der can overcome Anderson localization to facilitate ef- 
ficient transfer. However, we have shown that at higher 
hierarchy levels, dynamic disorder is reduced while static 
disorder is increased. Anderson localization therefore 
provides the ultimate limit to coherent long-range ex- 
citon transfer. Transfer over scales longer than dk A must 
be by thermally induced hopping. 



Conclusions 



over scales larger than achievable in equivalent unclus- 
tered networks. 

Anderson localization due to static disorder provides 
the ultimate limit on coherent long-range transfer, rather 
than thermal dynamic disorder or large separations be- 
tween chromophore clusters. This suggests that deter- 
mining the significance of coherence-assisted transfer over 
physiologically-relevant length scales in real systems may 
require a focus on the scaling of static disorder, for exam- 
ple by studying inhomogeneous broadening. Beyond bio- 
logical systems, our results suggest the possibility of long- 
range coherence-assisted transfer at ambient tempera- 
tures in a chromophore network engineered with small 
static disorder and clustering geometry which minimizes 
decoherence at larger length scales. This may also hint 
at a more general principle concerning networks of elec- 
tronically coupled quantum systems in condensed matter: 
Hierarchically clustered network geometries, carefully de- 
signed to complement the dielectric and mechanical prop- 
erties of their surrounding media, may aid in sustaining 
long-range coherent interactions at high temperatures. 
This is potentially of interest for quantum information 
applications. 

Avenues for extension of our model include generaliz- 
ing the dynamical treatment to strong system-bath cou- 
pling at low hierarchy levels; reference 15j and refer- 
ences therein describe such non-Markovian treatments. 
We speculate that including dissipative processes and 
analyzing their scaling may be useful in studying non- 
photochemical exciton quenching mechanisms, in which 
long-range LHC clustering geometry is implicated [l7j . 
Finally, it may be of interest to generalize the renor- 
malization procedure to geometries which are higher- 
dimensional, have non- uniform site separations and/or 
have variations in cluster size between hierarchy levels, 
and to accommodate nonlinear phonon dispersion due to 
material differences between hierarchy levels. 



We have analyzed scaling trends in exciton transfer 
in a chromophore network with self-similar, hierarchi- 
cal, dimeric clustering geometry, and with Markovian 
system-bath interaction. We have assessed bath-induced 
crossover from coherent to incoherent transfer across the 
hierarchy's spectrum of length scales via a novel reno- 
malization procedure. For networks without static dis- 
order in the site energies, the transfer coupling varies 
inversely with both hierarchy level and degree of chro- 
mophore clustering. Counterintuitively, excitonic deco- 
herence also varies inversely with both level and cluster- 
ing. In a biologically relevant parameter space, our cho- 
sen network geometry supports coherent transfer for lev- 
els 0-3, and also higher levels in more clustered networks. 
Chromophore clustering therefore facilitates long-range 
coherence rather than inhibiting it. This may suggest 
that the clustering hierarchies ubiquitous in biological 
chromophore networks help exploit excitonic coherence 



SUPPLEMENTARY INFORMATION 
Model 

We assume for the complete Hamiltonian (20j ] . 

H = H s + H B + H SB , where (1) 
H s = i?( 1 )|l)(l| + i?( 2 )|2)(2|-A(|l)(2| + |2)(l|),(2) 

H B = ft]>>(qMq)ta(q), (3) 
q 

h sb = ? ;^G( 12 )( q )(a( q )- a (- q )t). (4) 
q 

Here, {|1), |2)} is the single-excitation basis for two sites 
with energies E^' 2 \ static disorder e = E^ 1 ' — E^, and 
transfer coupling A. a(q) annihilates a phonon with wave 
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vector q, G^ 12 ^(q) = M(q)p(q) is an effective exciton- 
phonon coupling operator for the dimer. 



System-bath coupling and spectral density 

The spectral density depends on the microscopic 
system-bath coupling, g^ 12 \q). Details of system-bath 
couplings for biological LHCs are not yet well known 
11 1 , but are nonetheless generically described [20] by an 
effective exciton-phonon coupling operator, & 12 \q) = 
Af(q)/5(q), with M(q) the coupling strength and 



p(q) = W/d 3 r^V)^(r) 

m,n \J 

= P(q)^e- tq - r '" \m){m 



<1 *^ I C^Cn 



(5) 



the Fourier transform of the exciton density operator. 
ipjr(r) is the wave function for an excitation on site i, 
and 'P(q) = J (i 3 r|-0x(r)| 2 e~ 4q ' r is the form factor of 
the wave function [24( . We assume single-site excitations 
are localized within length I much smaller than the in- 
tradimer separation d, such that 'P(q) ~ 1 for all |q| 
values of interest. Evaluating equation <j5j) for the dimer 
comprising sites 1 and 2 gives, 



p(<D 



cos(^) (|1)(1| + |2)(2|) 

+ isin ^_ff) (|1>(1|-|2)(2|) 

'd • q\ 
sm | — - — ] Z, 



where Z — |1)(1| — |2)(2|. The term proportional to the 
identity is ignored since it only shifts the phonon energies, 
and the phase factor is also ignored. Then we can write 
G*(i2)( q ) =5 (i2)( q ) z, where 



ff ( 12 )(q) = A/(q)sin 



(6) 



The system-bath interaction Hamiltonian is therefore 



H SB = ZY,9 (12) ^) (a(q)-o(-q)t 
q 



(7) 



In order to compute the spectral density, we assume as 
a first approximation, deformation coupling of the form 
M(q) = D\q\y l /h/2/j,Vw(q), which accurately describes 
a quantum dot in a crystal [2(J. Here, p is the mass 
density, D the deformation potential, and V the quan- 
tization volume. We stress that the precise form of this 
coupling is used as a guide only, and the salient features 
of our model on which we base our conclusions are inde- 
pendent of the specific parameters therein. The spectral 



density at the Oth level is computed by assuming linear 
phonon dispersion in the bath and approximating with 
the integral, 



V 



J{uj) = — I dqd(t>d6q*sm<p5(u-vq)\g ( ™] g \ 2 , (8) 



4tt 2 



which, using equation evaluates to equation (1) in 
the main text. 



Dimer master equation and dynamics 

Here we derive a non-dissipative Markovian master 
equation for pi, the combined density matrix of the 
dimer site populations in an interaction picture with 
respect to H s + H B . Non-Markovian effects enhance 
the preservation of excitonic coherence at intra-LHC 
scales, compared with the predictions of Markovian mod- 
els (6 " 



Hi lla [25( • However, the ultimate focus of this 



analysis is long-range transfer, which we show to be less 
affected by system-bath interactions than short-range 
transfer, and therefore more appropriate for a Markovian 
description. Moreover, a Markovian treatment at intra- 
LHC scales is a conservative choice since it predicts faster 
decoherence than non-Markovian treatments [3 [Tj| . 

Transforming equation (J7J into an interaction picture 
with respect to H S +H B , in which a(q, t) = a(q)e~ iLJ ( q ) t , 
we find 



P{uj')e- luJ,t + P(u')U iu ' t , 



V4A 2 



The 



operators 
and 



)o+ + (e + 0)ct_ + 2AZ) , where 



where <j> = 

m = % 
m = ^(( 

a+ = |1)(2|, a- = |2)(1|, Z = |1)(1| - |2)(2|. The terms 
proportional to Z cause pure dephasing in the site basis 
while those proportional to <r + or <r_ induce transitions 
between sites. The term proportional to the identity can 
be ignored, giving the interaction picture Hamiltonian, 

#/ SB (t) = Z(t)E<? (12) (q) (a(q)e~ Mq)t +a(q)V" (q) * 
q 

Integrating over the von Neumann equation for the den- 
sity matrix, W, of the total system then tracing out the 
phonon modes gives for the combined density matrix of 
the site populations, 



Pl (t) = / dt'Tr ph {[Hf B (t), [Hf B (t'),Wj(t')}}}. 

J to 

We assume weak system-bath coupling and rapid bath 
relaxation, allowing a Born-Markov approximation to be 
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made by setting to 
The result is 



-oo and replacing Wi it') — > Wiit), with eigenenergies £ 



(ij) 
k± 



). Here, 



pi= W)wv) + vd[pV)\pi 

u)'=0,<f> 

+N{uj')V[P{uj , i i ]p I ), 

where V[A]p = ApA^ - \{A^ Ap + pA^A). N(uS) = 
^ e hui/k B T _ ^-1 j g ^ e ph onorL occupation of the bath 
at temperature T and J(ui) is the spectral density. Ex- 
citon recombination processes are not included, since we 
assume the lifetime to be long compared with the time 
required for transfer to a RC, which is necessary for the 
near-unity quantum efficiencies of unstressed photosyn- 
thetic membranes. 

In the ordered case (e = 0), an initially localized exci- 
tation (pi i 1 (0) = 1) evolves in the Schrodinger picture 
as 



pSiAt) 



lie 



-F(2A)t 



cos(2At) 



J 



1,2, 



showing coherent tunnelling between sites at frequency 
2 A, damped by bath-induced decoherence at rate 
F(2A) = J(2A)(JV(2A) + i), where N(lu) = ( e hu / kBT - 
is the phonon occupation of the bath at temperature 

T. 



Renormalization procedure 

Renormalization from the fcth to the (fc + l)th hierar- 
chy level is achieved via a coarse-graining operator which 
maps the n single-excitation electronic eigenstates (exci- 
tons) of a chromophore cluster at the fcth level to the 
unique single-excitation eigenstate of the equivalent ef- 
fective site at the (k + l)th level. To construct the oper- 
ator, we consider a pair of dimers (tetramer) (Figure 2 in 
the main text) for which the Hamiltonians ([I]), ([2]) and 
(|4j) are generalized as 



H, 



(1-4) 



H, 



s(l-4) 



SS(!- 4 ) 



i i,j 

for i ^ j; i,j = 1..4, 
'E 6 o J> (q) (o(q)-o(-q) f 



H, 



with numerical subscripts labeling hierarchy level. The 
bath Hamiltonian H B is still given by equation ([3]) . For a 
dimer comprising sites i and j at the fcth hierarchy level, 
with effective static disorder e 
and transfer coupling A 

|±)«fl = co S e k ^\i) k + S m^\j)k}, where k ™ 
arccos (-(e™ ± ^)/^A k ^ + (J) ± ^))2 



= E 



(i) 



E 



(j) 



(ij) 



the eigenstates are 



0®) = \jA/S^ iY + e^' . The interdimer transfer 
dynamics are dominated by the low-energy eigenstates 
of the tetramer. Therefore, wc define a renormalization 
operator for the dimer, which projects onto its lowest 
energy eigenstate. For j, n even and i = j— 1, m = n— 1, 



MY 



^ = li/2) 



(9) 

The complete Hamiltonian for a dimer at the (k + l)th 
level is obtained by renormalizing the corresponding 
tetramer at the fcth level: 



H, 



(12) 
fc+1 



^^(34)^.-4) ^^(34)1 



= H, 



S (12) 



H k+1 



(12) 



E. 



fc+1 

(1) 
fc+1 



SB (I 2 ) 
fc+1 ' 



where 



l) fe+1 (l|+4 2 4 ) 1 |2) fc+1 (2| 



Ai 1 + 2 l(|l) fc+1 (2| + |2) fc+1 (l|), 



q 

(«7) 

-ffs commutes with pj, , all levels of the hierarchy are 
coupled to the same bath, although the coupling, de- 
scribed by equation (3) in the main text, changes with k. 
Since the complete Hamiltonians at the fcth and (fc + l)th 
levels have the same form, the dynamics at both levels 
are described in the same way, with only the effective 
site, transition coupling and system-bath coupling ener- 
gies being transformed between levels. 



Renormalization of system-bath coupling and 
spectral density 

The spectral density is not renormalized directly but 
rather, at the fcth hierarchy level it is computed from 

( 12 ) 

a renormalized system-bath coupling g\ (q). Equation 
([5]) is renormalized using the operator defined in equation 
©, letting the subscript — > k for generality: 

n (r,\ - R [1 ^R^n ( n\ «(l 2 ) t «(34) t 

Pfc+uqj-Pfc Pfc Pfc(qj pfc Pfc 

COS 



((exp (id k+1 • q) cos 2 9^ + cos 2 0f_ 4) ) I fe+1 
exp (id k+1 ■ q) cos 2 0[__ 2) - cos 2 0f_ 4) ) Z k+1 



cos 



dfc • q 



dfc+i • q 



,c2fl(!2) 



'k- 



e 2 fl (34) 
,S d k- 



( dfc+1 ' q\ / 2 fl (12) 

sin I I [cos 8 k _> 



cos 2 e 4) ) 



Zk+l- 



Therefore, the system-bath coupling for the renormalized 
dimer at the (fc + l)th level is given by equation (TT5)) . For 



7 



1 .0 




/ \ 1 




0.8 








0.6 








0.4 




vV \ 2 




0.2 




fl\ 3 




0.0 









5 10 15 

oj QO'V) 



FIG. 1: Spectral density modulating functions, shown 
for first five renormalized hierarchy levels k (curves labelled). 
Clustering is set at di/do = 3.4, in analogy to trimeric LHCII. 



At higher levels of the hierarchy, the spectral density 
(equation (4) in the main text) contains a modulating 
function which depends on the network geometry: 

1 f 1 

W k(u) = 2k+T J dx[(l - ax(A k (u)x)) 

k-l 

x JJ(l + cos(A 3 -(w)x))]. (10) 

3=0 

Figure [T] plots Wfc(w), for the first five renormalized lev- 
els, assuming d\/do = 3.4, in analogy to trimeric LHCII, 
the major LHC in higher plants, where the n-n intra- 
monomer chromophore separation do ~ 10A and the 
centre-centre intermonomer separation eJi w 34A. 
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FIG. 2: Variations in the decoherence rate, Fk(2Ak), 
with degree of network clustering (di/do), for first six hierar- 
chy levels k (curves labelled). Temperature, T=293K. 



an ordered network (eg = 0), equation (|T5|) reduces to 



<0(q) = M(q)cos 



d fc • q 



dfc+i • q 



and recursive renormalization gives equation (3) in the 
main text. The spectral density at the fcth (k > 1) hierar- 
chy level, Jfe(w), is then found by letting g^ 12 \q, </>, 9) — > 
9k (?j 4>i 6) m equation ((SJ) and substituting equation 
(3) from the main text. The result is equation (4) in the 
main text. In the small frequency limit where u) -C v/dk, 
with v the speed of sound, equation (5) in the main text 
allows for direct transformation between spectral densi- 
ties at different hierarchy levels. Outside this limit, how- 
ever, renormalization can proceed only via the system- 
bath coupling. This suggests that investigating multi- 
scale exciton transfer in a specific LHC network or other 
real chromophore network may require the determination 
not only of the spectral density at some scale, but also 
of the exciton-phonon coupling parameters at the scale 
of individual chromophores. 



Hierarchy energy transformations 

Through recursive renormalization we obtain exact ex- 
pressions for the effective site, transition coupling and 
system-bath coupling energies at the fcth (k > 1) hierar- 
chy level in terms of the 0th level energies in an ordered 
network (eg — 0): 



E 



(i) 



fc-i 



A 



(2m— 1,2 



(11) 



i=l 
2 k+i 



i=0 



E 



(2) 



. ipg) 



i=2 k +l 

k-l ^ 2 k ~ i 

^ / 2k-i-l ^ / 

i=0 m=2 k - i - 1 +l 

2 k p 2 k q 

l=2 fc O-l) + l j=2 k (q-l) + l 
k-l 



^(2m-L2m)^2j 



(ij) 




g™ (q) =Af(q)sin(^)n 

3=0 



dj-q 



(13) 
(14) 



Generalizing these transformations for a network with 
static disorder is nontrivial. Here we give transforma- 
tions between the £;th and (k + l)th hierarchy levels in 
a disordered network, which may be iterated to compute 
energies at an arbitrary level: 



E 



(i) 
fc+i 



E. 



(i) 



4A 



(12) J ,(12) 



4A 



(12)' 



(12) ,(12) 



k ) 2 



P (2) 
h k+l 



E 



(3) 



4A 



(34)" ,(34) 



4Ai 34) +(4 34) "4 34) ) 2 ' 
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A 



(12) 
fc+1 



(12) 



( 7 f )Af» > + 2A< 12 >A< 14 >) + ( 2 A< 12 >A 



ffiS(q) = M(q)cos 



(4A fc 
dfc • q 



(12)' 



• 7 f) 2 )(4Af) 2 



(34)' 
7fe 



(13) 

i/T" 



7 (3 4 ) A (23) 



COS 



dfc+i • q 



cos 



2/)(12) 



— COS 



e 4) ) 



sin 



dfe+i • q 



COS 



2/)(12) 



COS 



2 C 4) (> 



5) 



where 7^ 



Parameterization and scale-based crossover 

Parameterizations of the transfer coupling A range 
from more exact methods I26|. |27l to the simpler ideal 
dipole approximation (IDA) [26j. We implement the lat- 
ter as A® = K|7 ? W||77^|/(47re^ ) rf( ij ') 3 ), where 77^'^ 
are the transition dipole moments of sites i and j, and 
Cj, the dielectric permittivity between them, k is a fac- 
tor determined by the relative orientations of ^W'W and 
their connection vector. We set n = 1 for convenience, 
since our model accommodates arbitrary configurations. 
In real systems, jy( m ) varies if multiple chromophore types 
are present, and may change with each coupling, de- 
pending on the multiscale geometry of the network and 
material composition of its surroundings (eg. protein em- 



bedded in lipid, suspended in aqueous fluid) [16|,l28|. We 
avoid an elaborate, arbitrary choice of composition and 
geometry, assuming a single chromophore type with tran- 
sition dipole moment magnitude 77 = 4.51? [29[, and also 
uniform dielectric permittivity of ex — 2.1eo, where eo is 
the vacuum permittivity, as for intracomplex couplings 
in LHCII :'>() . We take n = 1500kg /m 3 for the mass 
density and v = 2500ms -1 for the speed of sound. d is 
set at lOA, based on the n-n spacings in LHCII. To our 
knowledge, the deformation potentials for chromophores 
in biological LHCs are not yet known. We estimate the 
deformation potential to be of order D rj O.leV" since 
this number gives physically plausible decoherence rates 
(Figure and predicts transfer dynamics which fall nar- 
rowly inside the coherent regime at the 0th hierarchy level 
(Figure 3 in the main text), which is consistent with ex- 
perimental findings for LHCs [H-Q- 

An intuitive explanation for the scaling of the trans- 
fer dynamics is as follows. At the fcth heirarchy level, 
the network geometry and surrounding dielectric set the 
transfer coupling rate 2Afc, which determines the value of 
the modulating function W k (2A k ) by which the spectral 
density J k (2A k ), and accordingly the decoherence rate 
F k (2A k ), are modulated. The ratio of transfer rate to 
decoherence rate then determines whether the dynamics 
are coherent or incoherent. Unsurprisingly, A k varies in- 
versely with hierarchy level and degree of network cluster- 
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FIG. 3: Variations in the effective intradimer transfer 
couplings, 2Afc, with degree of network clustering (di/do), 
for first six hierarchy levels k (curves labelled). 



ing (di/do), though increased clustering curtails the rel- 
ative reductions in A k between levels (Figure EJ. Coun- 
tcrintuitively however, F k (2A k ) also varies inversely with 
level and clustering, with (for our chosen parameters) a 
large reduction between the 0th and 1st levels and smaller 
reductions thereafter (Figure [2]). For large fc, 2A k <C 
v/d k and F k {to) k, Bk B Td k uj 4 /(6v 2 ). Hence for our 
chosen geometry, F k+1 (2A k+1 ) w 4(A fe+ i/A fe ) 4 F /£ (2A /£ ). 
Therefore, for sufficiently large k, the decoherence rate 
vanishes. 
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